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SUMMARY 


This  report  describes  the  research  performed  during  two 
years  of  support  by  the  AFOSR  under  grant  number  AFOSR-86-01 36 . 

The  research  falls  into  four  distinct  categories: 

(i)  A  continuation  of  the  Pi's  investigation  of  the  theory 
of  structure  functions  on  continua 

(ii)  The  construction  of  a  diffusion  model  for  a  system 
subject  to  continuous  wear 

(iii)  The  introduction  of  criteria  for  reliability  growth 

(iv)  An  investigation  of  the  stability  of  stochastic  models 
The  results  obtained  are  described  in  detail  in  Sections  1-4  below. 


□  □ 


CONTINUUM  STRUCTURE  FUNCTIONS 


A  continuum  structure  function  (CSF)  is  a  mapping 
y:A  ->[0,1]  ,  where  A  =  [0,1]”  ,  which  is  nondecreasing  in  each 
argument  and  which  satisfies  y(0,...,0)  =  0  and  y(l,...,l)  =  1. 

The  theory  of  such  functions  generalizes  traditional  binary  and 
multistate  structure  function  theory,  permitting  more  realistic 
and  flexible  modelling  of  systems  subject  to  reliability  growth, 
component  degradation  and  partial  availability.  CSFs  were 
introduced  by  Block  and  Savits  [9]  and  Baxter  [2]  and  subsequently 
studied  in  [3-7],  [13-15]. 

Define  P^  =  [x|y(x)  >  a  whereas  y(y)  <  a  for  all  y  <  x]  where 
y  <  x  means  that  y  <  x  but  that  y/x(0<a<l).  Block  and  Savits 
[9]  prove  that  if  y  is  right-continuous,  i.e.  if  each  Ua  =  {x|y(x)  > 
a}  is  closed,  then  each  Pa  is  nonempty  and  that  the  Pa ' s  characterize 
y  by  means  of  the  decomposition 

y(x)  =  max  min  I(x  >y  ,dcx 
yeP^  1 <i <n  i ■  i 


where  I  denotes  the  indicator  of  A. 

A 

Structure  functions  for  which  each  P  is  finite  are  of 

a 

particular  interest  as  we  shall  show.  Intuitively,  in  two-dimen¬ 
sional  space,  if  Pa  is  finite,  each  segment  of  3Uo  should  be  parallel 
to  one  of  the  axes.  This  is,  in  fact,  the  case  in  n-dimensional 
space . 
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NOTATION  Bc(x)  denotes  an  open  ball  of  radius  e  centered  at  x 


U(x)  denotes  {y|y  >  x} 

Theorem  (Characterization  of  Finite  P  ) 

—  •  ■*■  ■  oc 

Suppose  that  Ub  is  closed  (0  <  a  <  1).  Then  Pb  is  finite  if  and 
only  if  there  exists  an  e  >  0  such  that  Bs(x)nUQ  =  BE(x)nu(x)  for  all 
xeP  . 

a 


The  class  of  right-continuous  CSFs  for  which  each  P^  is  finite 
may  be  characterized  as  a  natural  generalization  of  the  class  of 
Natvig  CSFs  [3],  the  defining  binary  structure  functions  of  the 
latter  being  replaced  by  multistate  structure  functions.  In  order  to 
define  this  class  of  structure  functions,  it  is  necessary  to 
introduce  some  terminology. 

Definition 

Let  Sx  =  {0,1,..., M. }  for  some  non-negative  integer  M. , 

1  n 

i  =  1,2, ...,n.  A  function  h:A->X  S.  is  said  to  be  a  reduction  mapping 

i  =  1  1  " 

if  h(x)  =  (h  (xx ) , . . . ,hn (xn )  )  for  all  xeA  where 

h1:{0,l]-»Si  is  a  surjective,  nondecreasing,  r  ight-continu  jus  step 
function  for  i  =  1,2, ...,n. 

For  each  zeS. ,  define  h~ 1 ( z )  =  {xe[ 0 , 1 ] | h. ( x )  =  z};  clearly, 

inf  h”  1  (  z  )  exists  and  h. (inf  h^1 (z) )  =  z.  It  can  be  shown  that 

inf  h~ 1 ( z )  exists,  is  unique  and  equals  (inf  h‘ 1 ( zy  ),..., inf  If  (*n>>- 

n 

A  mapping  X  S.  -»{ 0 , 1 .  .  .  ,  M }  ,  M  >  1,  which  is 
i  =  l  1 
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nondecreasing  in  each  argument  is  said  to  be  a  multistate  structure 


function  (MSP).  An  MSF  is  said  to  be  coherent  if  for  any  ieC  and 

n 

j  >  1,  there  exists  a  ze  X  S.  such  that  <f>((j-l).  ,z)  <  <f>(j.  ,z)  where 

i  =  l 

( j i #  z  )  denotes  ( zx , . . . , z . _ x , j , z . + x ,  .  .  .  ,  zn )  [12]. 

n  n 

For  a  reduction  mapping  h:A  -»  X  S .  ,  we  write  S.  =  X  S.  , 

i=l  1  i=l  1 

the  range  of  h,  and  Lh  =  {inf  h~ 1 ( z ) | z eSh } .  We  also  write 

xA  =  {x.ex|ieA},  hA ( x )  =  { h. ( x . ) ch ( x ) | i eA)  and  SA  =  {zA|zeSh} 

where  AcC . 


Definition 

Ca 

Let  {(♦a»h<sa),  0  <  a  <  1)  be  a  collection  of  pairs  of  functions 
such  that 

(i)  ha:£>  _>Sh  is  a  reduction  mapping  for  each  ae(0,l] 

c  « 

(ii)  ♦a:Sh->{0,l}  is  a  coherent  MSF,  the  components  of  which  are  the 

a 

elements  of  Co  for  each  ae(0,l] 

(  i i i )  U  Ca  ={l,2,...,n} 

a  C  C 

(iv)  For  each  xe  A,  <p  (h  “  ( x )  )  >  4>a  (h„e  (x)  )  whenever  a  <  3 

(0  <  a, 3  <  1 )  . 

The  function  is  said  to  be  an  F-type  CSF  if  it 

satisfies  the  condition 

C 

r  ( x )  >  a  iff  *  (h  a  ( x )  )  =1  (0  <  a  <  1) 

**  o  a 

for  all  xeA. 
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Definition  [ 3 ] 

A  CSF  y  is  said  to  be  weakly  cohe  rent  if 

sup  [rdj  ,x)  -  y(  0.  ,x)  ]  >  0. 
xcA 

Theorem 

A  weakly  coherent  CSF  is  of  F-type  if  and  only  if  it  is 
right-continuous  and  each  Po  is  finite. 

The  importance  of  F-type  CSFs  is  seen  in  the  case  where  X1  ,  .  .  .  ,  X  , 
the  states  of  the  n  components,  are  independent  random  variables. 

A  computationally  tractable  expression  for  <f>(x)  =  P{y(X)  >  x}  ,  the 
stochastic  performance  function,  occurs  only  in  certain  special  cases. 
If,  however,  Px  is  finite,  it  can  be  shown  that 


N  n  -  ,  -j  .  n_  U  .  >  (  j .  > 

♦  (x)  =  Z  Z  F.  (y.  31  )  -  Z  Z  H  F.(y.  1  Vy.  2  ) 

j"l  i  =  1  j  j  <  j  2  i-1 

n 

+  ...  +  (-1)N  1  n  F. (max  y.*'’1) 

i  =  l  1  l<j<N  1 

where  F .  ( t )  •=  P{X.  >  t }  ,  i  =  1,2,  ...,n,  so  that  <Hx)  can  be  easily 
evaluated . 

More  is  true.  If  y  is  an  arbitrary  right-continuous  CSF  with 
stochastic  performance  function  <f>,  there  exists  a  sequence  {  y„  }  of  CSFs 
such  that  pointwise  as  m-»®  where  $  is  the  stochastic  performance 

function  corresponding  to  y  .  Further,  if  <f>  is  continuous,  the 
convergence  is  uniform. 
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A  condition  under  which  <f>  is  continuous  is  given  by  the  following 


proposition . 

Definition 

Write  x  <  y  if  x.  <  y  and  xi  <  y.  if  y.  >  0,  i  =  1,2, ...,n. 

A  CSF  y  is  said  to  be  strictly  increasing  if  y(x)  <  y(y)  whenever 
x<y . 

Proposition 

If  r  is  a  strictly  increasing  CSF  and  if  X  is  an  absolutely 
continuous  random  vector  then  $  is  continuous. 

Suppose,  now,  that  the  structure  function  changes  with  time, 
possibly  reflecting  improvements  as  the  system  is  developed  or  changes 
in  the  use  of  the  system.  Let  rt  denote  the  CSF  at  time  t.  Suppose, 
further,  that  the  states  of  the  components  comprise  a  stochastic 
process  (X(t),  t>0],  perhaps  reflecting  the  degradation  of  the 
components  with  use  or  the  replacement  of  failed  components. 


Theorem 

If  X(t)^X,  an  absolutely  continuous  random  vector,  as  t-»® 
and  rt  -»r  pointwise  on  A  as  t-*®,  then  yt  (X(  t )  )^y(X)  as  t-»®. 

The  proof  is  by  consideration  of  the  set  E  =  {x  e  A | y  (xfc )^y( x ) 
for  some  class  { x t }  such  that  lim  xt  =  x} .  It  can  be  shown  that  if 

t-»oo 

Yt  *+y  pointwise  as  t->®,  then  E  has  Lebesgue  measure  zero  and  hence,  if  X 
is  absolutely  continuous,  PoX_1(E)  =  0  so  that,  by  Theorem  5.5  of 
Billingsley  [8,  p.  34],  the  result  follows. 
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Consider,  now,  the  class  {  Yt  }  of  CSFs  where  r't  ■  [  0 , 1  ]  -»[  0 , 1  ]  is 
defined  by 


fl  if  x  >  1/2  +l/t 

YJX)  - 

VO  otherwise 

for  t>2  and  the  CSF  y' :  [  0 , 1  ]-»(  0 , 1  ]  defined  by 

fl  if  x  >  1/2 

Y'(X)  = 

VO  otherwise. 

Clearly,  y^  does  not  converge  pointwise  to  y' .  Since,  however,  the  set 
of  points  (viz.  {1/2})  at  which  convergence  fails  has  Lebesgue  measure 
zero,  and  hence  PoX~  1  measure  zero  if  X  is  absolutely  continuous,  it  is 
reasonable  to  enquire  whether  there  exists  an  alternative  mode  of 
convergence  under  which  y^  (  X(  t )  )  ^y'  (  X )  as  t-»®.  This  motivates  the 
following  definition. 

Let  A  be  the  group  (under  composition)  of  all  homeomorphi sms 
X:  Rn  -»Rn  such  that  X(0,...,0)  =  0  and  A(l,...,l)  =  1  and  let  d  be  the 
Euclidean  metric  on  Rn .  For  any  two  right-continuous  CSFs  y  and  y2 
we  define 

S  ( y.  ,y,  )  =  inf  |  sup  |  y  (y*)  -  y,(X(y)*)|  +  sup  d(y*,X(y)*)i 
XeA  V  yeRn  yeRn  J 

where  y*  is  the  vector  whose  ith  component  is  (y. V0)A1. 

It  is  easily  verified  that  S  is  a  metric  on  the  space  of 
right-continuous  CSFs;  we  call  it  the  quasi-Skorohod  metric  and  we  call 
the  induced  topology  the  quasi-Skorohod  topology  (QST). 
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Convergence  in  the  QST  is  weaker  than  uniform  convergence;  it  is 
neither  weaker  nor  stronger  than  pointwise  convergence. 

Theorem 

Suppose  that  X(t)^X,  an  absolutely  continuous  random  vector,  as 
t -*<*>.  Suppose,  further  that  y  is  a  right-continuous  CSF  and  that  { y  } 
is  a  class  of  right-continuous  CSFs  such  that  y  ->y  in  the  QST  as  t -»«*>. 
Then  yt  (X(  t)  )%(X)  as  t->®. 

For  the  CSFs  { y' }  and  y'  defined  on  page  7,  it  can  be  shown  that 
yt'  (  X  ( t )  )  ^y '  (  X )  as  t-»°°  whenever  X(t)^X,  an  absolutely  continuous  random 
vector,  as  t->°°. 

Let  i>a  denote  the  intersection  of  8Uh  and  {cx|0  <  a  <  1},  the 
diagonal  of  A.  We  say  that  is  the  key  vector  of  Ua  and  we  call  the 
scalar  5a  the  cor responding  key  element .  Kim  and  Baxter  [13]  show  that 
for  any  CSF  y,  the  key  vector  always  exists,  is  unique  and,  if  y  is 

continuous,  y( )  =  a  for  all  ae(0,l].  Kim  and  Baxter  [13]  use  the  key 

element  to  define  reliability  importance  when  X  is  a  random  vector: 

they  define  the  reliability  importance  of  component  i  at  level  ae(0,l] 
as 

R.(a)  =  P  { y  ( X )  >  a  |  X.  >  SJ  -  P[y(X)  >  a|Xi  <  §a  }  , 

i  *=  1,2, ...,n.  This  generalizes  Birnbaum's  definition  of  reliability 
importance  for  the  components  of  binary  structure  functions  [1, 

Chapter  2]  to  the  continuum  case. 
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For  any  component  i  and  any  subset  Aca,  define  A1  =  {xeA|(-.,x)  = 
(•^z)  for  some  zeA}.  Notice  that  aca1  and  that  A  =  A1  if  and  only  if 
whether  or  not  xe A  does  not  depend  on  the  state  of  component  i.  Let  /v 
denote  Lebesgue  measure  on  Rn . 

Definition 

Let  y  be  a  CSF  and  1^  xe[0,l].  We  say  that  component  i  is 

almost  irrelevant  to  y  at  level  a  if  there  exists  a  subset  E  ca  such 

-  -  a 

that  u(  Ec  )  =  0  and  U  flE  =  (U  DE  HE  .  Further,  if  component  i  is 

a  a  a  a  a  oc  c 

almost  irrelevant  to  y  at  level  a  for  all  ae[0,l],  we  say  that  it  is 
almost  i r relevant  to  y. 

It  can  be  shown  that  component  i  is  almost  irrelevant  to  y  if  and 

only  if  there  exists  a  CSF  y'  such  that  y’  -  y  a . e . [ p ]  and 

sup [ y( 1, ,X )  -  y( o  ,x ) ]  =  0. 
x  s  A 

Before  proceeding  to  a  study  of  properties  of  Ri(a),  it  is 

convenient  to  deduce  sufficient  conditions  under  which 

0  <  P{X.  >  6  }  <  1  for  0  <  a  <  1.  It  can  be  shown  that  the  conditions 

that  y  is  continuous  at  (0,...,0)  and  (1,...,1)  and  that  the  support  of 

each  X.  is  the  unit  interval  are  sufficient  to  ensure  that 
1 

0  <  P{X.  >  S  }  <  1  for  0  <  a  <  1. 

Theorem 

Suppose  that  the  CSF  y  is  continuous  at  (0,...,0)  and  (1,...,1) 
and  that  X  , ...,Xn  are  independent,  absolutely  continuous  random 
variables,  the  support  of  each  of  which  is  the  unit  interval.  Then  for 
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all  «e( 0,1 )  and  each  i  =  1,2, R.  ( a)  =  0  if  and  only  if  component 
i  is  almost  irrelevant  to  y  at  level  a. 

Definition 

A  CSF  y  is  said  to  be  strongly  increasing  if  y(x)  >  Y(y)  whenever 
x.  >  y  for  i  =  1,2, ...,n. 

Define  a  function  6  :  [  0 , 1  ]  ->[  0 , 1  ]  by  6(a)  =  6  where  6  is  the  key 
element  of  Ua  .  We  call  6  the  key  function  of  y.  It  is  easily  seen 
that  6  is  well  defined,  nondecreasing  and  left-continuous,  and  is 
continuous  if  y  is  strongly  increasing. 

Theorem 

Let  y  be  a  strongly  increasing  CSF  which  is  continuous  at 
(0,...,0)  and  (1,...,1)  and  suppose  that  X1,...,Xn  are  independent, 
absolutely  continuous  random  variables,  the  support  of  each  of  which  is 
the  unit  interval,  then  R.(a)  is  continuous  cn  (0,1)  for  i  =  l,2,...,n. 

It  is  of  interest  to  determine  when  one  component  is  uniformly 
more  important  than  another,  i.e.  when  R.  (a)  >  R^  (a)  for  all  a,  c.f. 
Theorem  2.1  of  Natvig  [17]. 

Definition 

Let  y  be  a  CSF  and  let  6a  denote  the  key  element  of  Ua .  We  say 
that  component  i  is  connected  _in  series  ( parallel )  to  the  remainder  of 
the  components  if,  for  all  ae(0,l],  xeUa=>  (<=)x.  >  6o  . 
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In  the  special  case  when  y  is  a  binary  structure  function,  this 
definition  reduces  to  the  usual  definitions  of  being  connected  in 
series  and  in  parallel. 

Theorem 

Let  y  be  a  CSF  which  is  continuous  at  (0,...,0)  and  (1,...,1)  and 
suppose  that  Xx,...,X  are  random  variables  each  with  support  the  unit 
interval.  If  component  i  is  connected  in  series  to  the  remainder  of 
the  components  and  if  X.  <l  x  ,  then  Rz(a)  >  R..(a)  for  all  ae(0,l), 

j  *  i. 

Notice  that  it  was  not  necessary  to  assume  that  Xlf...,X  are 
mutually  independent. 

Theorem 

Let  y  be  a  CSF  which  is  continuous  at  (0,...,0)  and  (1,...,1)  and 
let  Xx,...,Xn  be  independent  random  variables,  each  with  support  the 
unit  interval.  If  component  i  is  connected  in  parallel  to  the  remainder 
of  the  components  and  if  X.  V  X.,  then  R.  (a)  >  R..  (  a)  for  all 

at ( 0 , 1 )  ,  j  t  i . 

The  algorithm  previously  introduced  to  calculate  <Ma)  can  easily 
be  modified  to  evaluate  R. ( a) . 
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2  .  A  DIFFUSION  MODEL 

This  section  introduces  a  model  for  a  system  whose  state  changes 
continuously  with  time,  perhaps  reflecting  continuous  wear.  It  is 
assumed  that  the  state  of  the  system  is  initially  x0>0  and  thereafter 
follows  Brownian  motion  with  negative  drift  and  an  absorbing  barrier  at 
the  origin.  It  is  further  assumed  that  the  state  of  the  system  is 
increased  by  a  repairman  who  arrives  according  to  a  Poisson  process  of 
rate  \>0;  if  the  state  of  the  system  when  the  repairman  arrives  exceeds 
a  threshold  a  (0<a<x0),  no  action  is  taken,  otherwise  the  repairman 
instantaneously  increases  the  state  by  a  random  amount  Y  where  Y>a 
almost  surely. 

Let  X ( t )  denote  the  state  of  the  system  at  time  t,  t>0,  let 
F ( x , t )  =  P[X(t)<x]  denote  the  distribution  function  of  X(t)  and  suppose 
that  the  corresponding  density  f(x,t)  =  3F(x,t)/9x  exists  for  x>0. 

Since  {X(t),  t > 0 }  is  a  diffusion  process  with  jump  discontinuities,  it 
is  a  special  case  of  the  model  analyzed  by  Feller  [11]  and  hence,  from 
equation  (27)  of  Feller  [11],  we  obtain  the  following  integro- 
differential  equation  (valid  for  x>0): 

(*)  y£-f(x,t)  =  jo2  — f  (x,  t)  -  U~  f(x,t)  -  Xl(xSot)  f(x,t) 

3  x 

+  X/*I{x_z<a) f(x-z,t)g(z)dz  +  AF( 0 , t ) g ( x ) 

where  //  and  a2  are  the  parameters  of  the  Brownian  motion  (/j<0  ,  o>0)  and 
g  is  the  density  of  Y. 
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Feller  [11]  has  obtained  the  general  solution  to  integro- 
differential  equations  of  the  form  (*),  but  his  solution  is  extremely 
difficult  to  apply,  and  hence  a  different  approach  was  adopted:  an 
explicit  formula  for  F(x,t),  0<x<a,  was  obtained  by  a  purely  probabi¬ 
listic  argument  which  does  not  make  use  of  (*),  and  these  results  and 
(*)  were  used  to  obtain  an  expression  for  the  Laplace  transform  of 
f(x,t).  The  distribution  function  of  X(t)  is 

F ( x , t )  =  JpB(x,t-u)e-x 1 t-u ’dH(u)  ( 0 <x <a ) 


where 


1-B ( x , t ) 


(  Z-iX-fjt  ) 2 

2a2  t  / 


L2±.a_AtJ.l|ld2  ( x>0  ) 

2  a2  t 


and  where 

<x> 

H( t }  =  W  (t)  +  I  W  *K<n) (t) 

X  ,  Ot  .a  ,  o 

o  n=l  o 

where  the  asterisk  denotes  Stieltjes  convolution  and  the  superscript 
denotes  n-fold  recursive  Stieltjes  convolution.  An  expression  for  K  is 
given  by  the  following  formula: 

K(t)  =  / «  j" W  {s)Xe'x ' t_s ’ dV(u)ds 

u  (X  u  ,  a 

where 
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V(x)  =  J"  J*  B(  (  x-y )  Aa,  t )  \e  xtg(y)dydt 


+  f*[B(x,t)  -  B(a,t) 3Xe_Xtdt 


and  where 


wa  b(t)  =  J*  — exp{- (  a-b+^ux  ) 2 /2ct2  x  }dx  . 
cti  2  jix3 

An  expression  for  F(x,t),  x>a,  cannot  be  obtained  by  an  analogous 
argument.  One  may,  however,  solve  (*)  to  yield  an  expression  for  the 
Laplace  transform  of  f(x,t).  This  is  uninformative  and  awkward  to  deal 
with  but,  for  the  stationary  distribution,  viz.  F(x)  =  lim  F(x,t),  an 

t-»» 

explicit  expression  may  be  obtained.  This  may  be  differentiated  to 
yield  (admittedly  clumsy)  formulae  for  the  stationary  mean  and 
variance . 

The  preceding  analysis  can  readily  be  extended  to  the  case  where 
the  repairman  arrives  according  to  a  stationary  renewal  process  in 
which  the  distribution  function  of  the  inter-renewal  times  is  Q,  say. 

Define  T0  =  inf{t|X(t)  =  Oj,  the  first  passage  time  to  state  0. 

We  deduce  expressions  for  the  distribution  function  of  T0  and  for  E(T0) 
and  we  derive  the  limiting  distribution  of  T0  as  <*-*«. 

Let  p  be  the  probability  that  when  a  repair  is  performed,  the 
system  is  at  state  0.  Clearly,  defining  P(t)  =  Q( u ) du/J" Q( u ) du , 

p  -  j;B(0,t)dP(t)/j;B(oc,t)dP(t) 

Observe  that  T0  satisfies  the  following  relation: 
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n  n  n 

x  E  s  +  £  Y.  =  S  +  £  S  with  probability  p(l-p) 

0  x0  i=l  1  x  o  i=l  i 

for  n  =  0,1,2,....  where  Y.  =  Y  for  all  i,  where  Sa  is  the  first 
passage  time  from  a>0  to  0  in  Brownian  motion  and  where  =  denotes 
equality  in  distribution.  It  is  readily  seen  that 

P{Sy  <  t}  =  /;Wy0(t)dG(y)  =  D  ( t )  ,  say, 

and  hence  the  distribution  function  of  T0 ,  L(t)  say,  is  given  by 

L ( t )  =  I  W  0 *D( n 1 ( t)p( l-p)n 
n=0  Xo' 

where  D(0,(t)  is  the  Heaviside  function  and  Wa  0  is  the  distribution 
function  of  Sa .  Now 

CD 

E ( T  )  =  l  [ E ( S  )  +  nE(Sy ) ]p(l-p)n  • 
n-0  0 

Since  E(Sa_b)  =  (b-a)///,  and  hence 
E(Sy  )  =  loE(Sy  ^dG^Y)  = 
it  is  easily  seen  that 

E(T0  )  =  -lpx0  +  m(  1-p)  ]/PfJ 
where  m  =  E( Y) . 
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Theorem 


As  a->®,  T0/E(T0)  converges  in  distribution  to  a  unit  exponential 
variate . 

Now  we  deduce  an  expression  for  n  (t  ,t  )  = 

P{X(t)  >  x  for  all  te[t1,t2]}. 

Consider,  firstly,  the  case  x  >  a.  The  first  passage  time  from 
X ( t,  )  to  x  is  equal  in  distribution  to  Sv ,  .  .  .  Hence,  if  W 

denotes  the  distribution  of  Sa_b,  the  first  passage  time  from  a  to  b 
in  Brownian  motion  with  parameters  p  and  a2  (a  >  b) 

MMM  -  CV.'MMVy-V- 

Consider,  now,  the  case  x  <  o.  Let  (Zx(t),  t  >  0)  denote 
Brownian  motion  with  parameters  p  and  a2  ,  an  absorbing  barrier  at 
x  and  initial  condition  Zx ( 0 )  =  a.  Let  Bx(y,t)  =  P{Zx(t)  <  y}  denote 
the  distribution  function  of  Zx(t)  and  let  T  x  be  the  first  passage 

time  from  y  to  x  in  the  process  {X(t),  t  >  0}.  It  can  be  shown  that 

the  distribution  function  of  T  ,  L  (t)  say,  is  given  by 

n«0 


where  p  =  /" B  ( x , t )dP ( t )/J" B  ( a, t )dP( t ) .  Hence, 


MW  -  I 
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It  is  henceforth  assumed  that  the  process  (X(t),  t>0}  is 
stationary,  and  that  the  repairman  arrives  according  to  a  Poisson 
process  of  rate  X.  Let  C:  denote  the  cost  per  visit  of  the 
repairman,  let  C2  denote  the  cost  of  a  repair  (this  may  depend  on  the 
distribution  of  Y)  and  let  C3  denote  the  cost  per  unit  time  of  the 
system  being  in  state  0.  We  calculate  C(X),  the  average  cost  per 
unit  time  over  an  infinite  horizon  for  a  given  arrival  rate  X.  To  do 
so,  it  is  convenient  to  define  the  following  random  variables: 

Nx  =  number  of  visits  by  the  repairman  in  a  cycle 
N2  =  number  of  repairs  performed  in  a  cycle 
T2  =  cumulative  time  spent  in  state  0  in  a  cycle 

where,  by  a  cycle,  we  mean  the  interval  between  two  successive 
instants  at  which  the  state  of  the  system  first  crosses  the  threshold 
a  from  above  following  a  visit  by  the  repairman  (recall  that  the 
sequence  of  such  instants  comprises  an  embedded  renewal  process). 

The  duration  of  a  generic  interval  is  denoted  T* .  It  is  clear  that 

N2  =  I(z(T)<a)  and  that  Ti  =  ^  o 1 1 z { t ) s  o } dt  where  T  is  an  exponential 
variate  with  parameter  X  and  (Z(t),  t > 0 }  denotes  Brownian  motion  with 
parameters  v  and  a2 ,  an  absorbing  barrier  at  the  origin  and  initial 
condition  Z(0)  =  a. 

It  is  easily  seen  that  E(NX  )  =  XE(T*).  Further,  E(N2)  = 

-  a  e  (  x  ) 

p { Z ( T )  <  a}  and  it  can  be  shown  that  E(T2)  =  e  /X  where 

6a  (X)  -  ( fj+  ( jj+2\a2  )1/2  ]/c2  . 

Let  c"(X)  denote  the  average  cost  per  cycle  for  a  given  arrival 
rate  X.  By  the  foregoing, 
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-O0  (XI 

C  {  X )  =  C:XE(T*)  +  C2q(X)  +  c3  e  1  /X. 


Then  C(X)  =  C(X)/E(T*)  and,  since  E ( T* ) 


it  follows  that 

C(  X)  =  Cx  X  + 


-  a  e  (  x  ) 

C2  Xq(  X)  +  C3  e  1 

-  a  e  (  X  ) 

e  -  Xmq (  X) //j 


where 


-  a  0  (  X  ) 

e  1  /X  -mq (X)/p, 


q  (  X ) 


e2  (X)  -e1  ( x ) exp { a[ ©2  (X)-e1  (X)  n 

e2  ( X )  -  077x1 


and  where  02(X)  =  [  p-  ( p+2X<x2  ) 1 7  2  ]/a2  ,  for  X>0;  by  continuity,  we 
define  C ( 0 )  =  C3 . 

We  seek  to  minimize  C(X)  by  varying  X. 

Firstly,  we  consider  the  case  C3  =  0. 


Theorem 

If  C2  >  -mC 3/fj,  then  C(X)  achieves  its  minimum  value,  C3  ,  at 
X=0,  otherwise  C(X)  decreases  monotonically  to  -jjC2/ m  as  X-»®. 

Now  we  consider  the  general  case. 


Theorem 

If  C3  +  C2  >  then  C(X)  achieves  its  minimum  value,  C3  , 

at  X-0,  otherwise  there  exists  a  unique  X*  (0<X*<®)  which  minimizes 
C(X)  . 
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For  the  case  a=0,  we  can  obtain  an  explicit  expression  for  X* . 


Theorem 

Suppose  a=0.  If  Cj  +  C2  >  -mC3///,  then  C(X)  achieves  its 
minimum  value,  C3  ,  at  X=0,  otherwise  C(X)  is  minimized  at 

C3  -  y-C1  (c2+m C3//i) 

^  “  mC3  /fj 
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3 .  RELIABILITY  GROWTH 

Let  { N ( t ) ,  t > 0 }  be  a  stochastic  point  process;  we  identify  the 

points  of  (N(t),  t > 0 }  with  the  failures  of  a  system.  Our  concern 

is  with  systems  which  exhibit  reliability  growth ,  i.e.  systems 

whose  reliability  increases  in  time.  We  model  such  systems  by 

stochastic  point  processes,  the  points  of  which  become  "scarcer" 

according  to  various  criteria  introduced  in  this  section.  Define 

the  intensity  function  h(t)  =  lim  P{N(t+S)  -  N(t)>0}&_1  and  say  that 

S->0  + 

{ N ( t ) ,  t > 0 }  is  orderly  if  P{N(t+S)  -  N(t)>l)  =  o( 5)  for  all  t>0.  Let 
H ( t )  =  E [ N ( t ) ]  denote  the  expected  number  of  failures  in  (0,t],  let 
X  denote  the  time  from  the  (n-l)th  failure  to  the  n*1  h  failure 

n 

(n>2),  X  being  the  time  to  the  first  failure,  and  define 

N  (  t )  + 1 

the  forward  recurrence  time  at  t  as  y(t)  =  E  X.  -  t. 

i-1 

The  following  four  criteria  for  reliability  growth  are  posited: 

Cl  N ( t+s )  -  N ( t )  is  stochastically  decreasing  in  t  for  all  s>0 
C2  y(t)  is  stochastically  increasing  in  t 
C3  Xn  is  stochastically  increasing  in  n 
C4  H(t)/t  is  nonincreasing  in  t. 

The  relationships  between  these  criteria  are  as  in  the  following 
diagram;  no  other  implications  hold. 


Cl 

4 

C2  C4 

C3 


*for  orderly  processes  only. 


By  way  of  illustration,  these  criteria  are  applied  to  three 
models:  the  initial  faults  model,  the  nonhomogeneous  Poisson  process 

(NHPP)  and  the  renewal  process. 

The  Initial  Faults  Model 

Suppose  that  there  are  initially  N  faults  in  the  system,  each  of 
which  will  eventually  give  rise  to  a  failure.  The  times  until  the 
failures  actually  occur  are  assumed  to  be  independent  random 
variables,  each  with  distribution  function  F,  so  that  N(t)  ~ 

Bin ( N , F ( t ) ) .  The  sequence  of  times  at  which  a  failure  occurs  is 
distributed  as  the  order  statistics  corresponding  to  N  independent 
observations  of  a  random  variable  with  distribution  function  F.  It 
can  be  shown  that  a  sufficient  condition  for  Cl,  C2  or  C3  to  hold  is 
that  F  is  DFR ,  and  it  is  obvious  that  C4  holds  if  F(t)/t  is 
nonincreasing  in  t. 

The  Nonhomogeneous  Poisson  Process 

Let  { N ( t ) ,  t>0}  be  an  NHPP  with  hazard  function  X(t)  = 
dE[N(t)]/dt.  A  sufficient  condition  for  Cl,  C2  or  C3  to  hold  is  that 
X  is  nonincreasing,  i.e.  that  F,  the  distribution  function  of  the 
time  to  the  first  failure,  is  DFR.  Further,  C4  holds  if  F  is  DFRA. 
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The  Renewal  Process 

Brown  [10]  proves  that  a  renewal  process  in  which  the 
distribution  of  the  inter-renewal  times  is  DFR  satisfies  Cl,  C2  and 
C 4.  Trivially,  any  renewal  process  satisfies  C3. 

Theorem 

Criteria  Cl,  C2  and  C4  are  preserved  by  the  formation  of  series 
systems  whereas  criterion  C3  is  not. 
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STABILITY  THEORY 


4  . 


No  stochastic  model  is  a  true  representation  of  an  actual 
physical  phenomenon:  assumptions  such  as  independence,  equality  in 
distribution,  exponentiality  etc.  are  seldom  satisfied  in  practice. 

It  is  therefore  of  interest  to  investigate  the  stabi 1 i ty  of  such 
models,  i.e.  to  determine  how  a  violation  of  the  assumptions  affects 
the  properties  of  the  model.  Let  2^  denote  the  class  of  all  "input" 
to  a  stochastic  model  and  let  'if  denote  the  corresponding  class  of 
"output".  Thus,  for  example,  in  a  GI/G/1  queue,  Ue  V.  would  comprise 
the  arrival  process  and  the  sequence  of  service  times  whereas  Vclf 
w' uld  be,  for  example,  the  sequence  of  waiting  times.  Letli*  be  the 
class  of  "perturbed"  "input"  and  letl^  *  be  the  corresponding  "output" 
class.  Then  a  stability  analysis  would  be  to  determine  conditions 
under  which 


p  (  U,  (J*  )  <e  =>  p(V,V*)  <  6(e) 

where  e>0,  S(e)>0  and  p  is  a  probability  metric,  typically  the 
uniform  metric,  on  the  space  of  distribution  functions  (for 
simplicity,  we  write  p(X,Y)  to  denote  the  distance  between  the 
distribution  functions  of  X  and  Y). 

We  now  present  two  stability  analyses:  one  of  a 
characterization  of  the  bivariate  Ka rshall-Olkin  distribution  and  one 
of  the  estimation  of  the  mean  of  the  (univariate)  exponential 
distribution . 

Recall  that  if  G(x,y)  =  P{X: >x,X2 >y}  denotes  the  bivariate 
survivor  function  of  (X  ,X2),  the  Marshall-Olkin  distribution  is 


defined  by 


G (  x  ,  y  )  =  expf  X2  x  -  X,  y  -  X,  2  max  (  x  ,  y  )  } 
for  X  ,X2>0,  X  2 >0  and  x,y>0  (e.g.  Marshall  and  Olkin  [16]). 

Let  B  denote  the  class  of  all  bivariate  survivor  functions  of 
pairs  of  nonnegative  random  variables.  For  GeB,  define  the  hazard 
vector  (h1(t),  h  ( t ) )  =  V  [ -logG( t , t ) ] ,  assuming  that  this  exists,  and 
write  H  (x,y)  =  - BlogG ( x , y ) /3x ,  H2(x,y)  =  -  3 logG ( x , y ) /3y  and  c  = 
H1(0,0)  +  H  (0,0)  .  The  corresponding  bivariate  Marshall-Olkin 
distribution  is  defined  as 

^  fexp{-cy  -  h1  (  0  )  (  x-y ) }  if  x>y 

G(x,y)  =  | 

lexp  { -cx  -  h2  (  0  )  (  y-x  )  }  if  x<y, 

i.e.  G  is  the  bivariate  survivor  function  of  the  Marshall-Olkin 
distribution  with  Xx  -  c  -  h2  ( 0 )  ,  X2  =  c  -  h2  ( 0 )  and  Xx  2  = 
hx ( 0 )  +  h2(0)  -  c.  Notice  that  if  (Xx,X2)  has  survivor  function  G, 
then  E(X.  )  =  1/h.  (0)  and  E(xJ)  =  2/[ h.  ( 0 ) ] 2  (i-1,2)  and  F(X2X2)  - 
[l/h1 (0)  +  l/h2 (0) ]/c. 


Definition 

Suppose  that  GeB,  the  survivor  function  of  (X2,X2),  satisfies 
the  inequalities 

(i)  G(x+t,y+t)  <  G(x,y)G(t,t)  for  all  x,y,t>0 

(ii)  G. (x+t)  <  Gi(x)G. (t)  for  all  x,t>0  where  G.  is  the  survivor 
function  of  XA  (i=l,2). 
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Then  G  is  said  to  be  weakly  bivariate  NBU  ( WBNBU ) .  Further,  G  is 
said  to  be  WBNBU*  if  G  is  WBNBU  and  if  the  partial  derivatives 
9G(x,y)/3x  and  3G(x,y)/3y  exist. 

Define 

(3(G)  =  { [  1/hj  (  0  )  -  E(Xt  )]  +  [  l/h2  (  0  )  -  E(X?)) 

+  (2E(X1)  -  h1(0)E(X12)]  +  [2E(X2  )-h2  (O)E(X2  )  ] 

+  4[E(X1)  +  E(X2)  -  cE(X1X2)}/c. 


Theorem 

If  G  is  WBNBU*,  then  G  is  Marshall-Olkin  if  and  only  if 
(3(G)  =  0. 


We  now  quantify  the  stability  of  this  characterization  by 
deriving  a  bound  on  the  uniform  distance  between  an  arbitrary  WBNBU* 
survivor  function  G  and  the  corresponding  Marshall-Olkin  survivor 
function  G. 

Theorem 

p(G,G)  <  (l+c2ec )(e(G) ]1/3  for  all  GeB. 

Suppose  that  Y  is  exponentially  distributed  with  probability 
density  function  (p.d.f.) 
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f¥  (X)  =  i  e"x/"  (  x  >  0  ;  >  0  )  . 


To  perform  hypothesis  tests  on  p ,  one  makes  use  of  the  fact  that,  if 

Y  , Y  , ...,Yn  are  n  independent,  identically  distributed  random 

n 

variables,  each  with  p.d.f.  f  ,  then  2  l  Y i/p  ~  X]  n  •  In  practice, 

i  =  1 

the  assumption  of  exponent i al i ty  is  only  an  approximation;  it  is 

therefore  of  interest  to  enquire  how  well  the  x\  distribution 

n 

approximates  that  of  2  I  X,///  where  Xi  ,X2  ,  .  .  .  ,Xn  are  n  independent, 

i  =  l 

absolutely  continuous,  identically  distributed  random  variables  with 
common  distribution  function  Fx  and  mean  p  and  Fx  e*?  ,  the  class  of 
all  possible  perturbations  of  the  exponential  distribution.  If 
is  the  union  of  the  class  of  HNBUE  distributions  and  the  class  of 

/M 

HNWUE  distributions,  a  simple  bound  on  the  uniform  distance  p a 
n  n 

between  Zn  *  2  £  Y  /p  and  Wn  =2  I  X. /p  can  be  obtained  without 
i=l  1  "  i=l  1 

making  any  assumptions  concerning  the  mechanism  generating  the 
perturbation.  It  can  be  shown  that 


p(W  ,Z  )  <  — - —  M2/3  | 1 -(a/p)2  |1/3 
^  n  '  n  '  —  2 1  /  3  n  7  1 


where  o  =  var(X)  and 


n  -  1  -  (  n  -  1  ) 


M  = 

n 


/n ( n-1 ) n  e 
i ( n-l ) ! 


Thus,  if  X  is  HNBUE  or  HNWUE  and  if  the  coefficient  of  variation 

n 

of  X  is  close  to  unity,  then  the  distribution  of  2  E  X.  /p  is 

i  =  l 

uniformly  close  to  that  of  the  n  distribution.  In  particular,  the 
sampling  distribution  of  the  usual  test  statistic  for  hypothesis 
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tests  concerning  the  scale  parameter  of  the  exponential  distribution 
is  robust  with  respect  to  moderate  departures  from  exponentiality 
provided  that  the  perturbed  random  variable  is  HNBUE  or  HNWUE. 

Suppose,  now,  that  no  assumptions  concerning  J"  are  made.  It  is 
necessary  to  make  an  assumption  concerning  the  mechanism  perturbing 
the  exponential  variate  Y.  Three  possible  cases  are  considered: 

( i )  Mixing 

Suppose  that  an  exponential  distribution  is  "contaminated"  with 
an  arbitrary  non-negative  random  variable  with  distribution  function 
H ,  i  e  . 


Fx ( t)  =  (l-e)e't/x  +  e  H( t) 
where  X  is  chosen  so  that  //  =  (l-e)X  +  to  and  where 

a  =  J”  H(t)dt.  It  is  assumed  that  e  >  0  is  small.  It  can  be  shown 
that  P(wn , zn )  <  0( t  1/3 ) . 

( i i )  Additive  Error 

Suppose,  now,  that  X  =  Yx  +  Ze  where  Yx  is  an  exponential 
variate  with  mean  X  and  ZE  is  an  arbitrary  random  variable;  X  is 
chosen  so  that  fj  -  X  +  E(Zc).  It  can  be  shown  that 

p(Wn,  zn  )  <  (1  +  1/S) [ 5b( S) ]1/( 1  +  4  ’  [2E(  |ze  |  )  ]4/( 1  +  S ’ 

where  b(S)  =  E ( | X | 2  +  6  )  +  E(|Yx|2  +  i)  (6  >  0)  is  assumed  to  be  finite. 
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( i i i )Right  Censoring 

Lastly,  suppose  that  X  =  min(Yx,N)  where  N  is  a  non-negative 
random  variable  independent  of  Yx  and  where  X  is  chosen  such  that 
E ( X )  =  /j  .  It  can  be  shown  that 

p(Wn>  Zn  )  <  |X2-V  |  +  eX2  +  2c 

where  e  =  P{N<hj  and  where  h  is  the  solution  to  the  equation 
P{N  <  h}  =  ye~n/r ( h+v) , 
writing  y  =  max  (X,/y). 

In  each  of  the  above  cases,  it  is  clear  that  a  small 

perturbation  of  Y  (i.e.  p(X,Y)  is  small)  yields  a  sampling 

n 

distribution  of  2  I  X.  //j  which  is  uniformly  close  to  that  of  the 

i  =  l 

distribution. 

2  n 
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